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ABSTRACT 

By using the method presented by Isobe et al. (2002), the non-dimensional recon- 
nection rate Vi n /vA has been determined for the impulsive phase of three two-ribbon 
flares, where Vi n is the velocity of the reconnection inflow and va is the Alfven ve- 
locity. The non-dimensional reconnection rate is important to make a constraint on 
the theoretical models of magnetic reconnection. In order to reduce the uncertainty of 
the reconnection rate, it is important to determine the energy release rate of the flares 
from observational data as accurately as possible. To this end, we have carried out one 
dimensional hydrodynamic simulations of a flare loop and synthesized the count rate 
detected by the soft X-ray telescope (SXT) aboard Yohkoh satellite. We found that 
the time derivative of the thermal energy contents in a flare arcade derived from SXT 
data is smaller than the real energy release rate by a factor of 0.3 - 0.8, depending 
on the loop length and the energy release rate. The result of simulation is presented 
in the paper and used to calculate the reconnection rate. We found that reconnection 
rate is 0.047 for the X2.3 flare on 2000 November 24, 0.015 for the M3.7 flare on 2000 
July 14, and 0.071 for the C8.9 flare on 2000 November 16. These values are similar 
to that derived from the direct observation of the reconnection inflow by Yokoyama et 
al. (2001), and consistent with the fast reconnection models such as that of Petschek 
(1964). 

Subject headings: Sun: flares — Sun: X-rays, gamma rays — Sun: magnetic fields 



1. INTRODUCTION 

Magnetic reconnection is now widely believed to be the mechanism of energy release in solar 
flares. The evidence of magnetic reconnection found by space observations includes the cusp-shaped 
post flare loop (Tsuneta et al. 1992), the loop top hard X-ray source (Masuda et al. 1994), the 
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reconnection inflow (Yokoyama et al. 2001), downflows above post flare loops (McKenzie & Hudson 
1999; Innes, McKenzie, & Wang 2003; Asai et al. 2004a), plasmoid ejections (Shibata et al. 1995; 
Ohyama & Shibata 1997, 1998), etc (see reviews by Shibata 1999; Martens 2003). Magnetic 
reconnection also plays an important role in various explosive phenomena in astrophysical, space, 
and laboratory plasmas (Biskamp 1992; Tajima & Shibata 1997; Priest & Forbes 2000). 

One important question in the physics of reconnection is what is and what determines the 
reconnection rate in a plasma with large magnetic Reynolds number such as solar corona. The 
reconnection rate is the amount of reconnected magnetic flux per unit time and given by %B, 
where V{ n is the velocity of reconnection inflow into the diffusion region and B is the magnetic field 
strength. It is more convenient and physically meaningful to define it in non-dimensional form as 
Ma = Vi n /vA, where va is the Alfven velocity. In the following we call this non-dimensional value 
the reconnection rate. It is important to determine the non-dimensional reconnection rate because 
theories of magnetic reconnection predict the non-dimensional value. 

In classical Sweet-Parker model (Parker 1957; Sweet 1958) the reconnection rate Ma is given 

by 

M A = R- 1/2 , (1) 

where R m is the Magnetic Reynolds number defined by the Alfven velocity (also called the Lundquist 
number) : 

Rm = (^), (2) 

r) 

L is the length of the reconnecting current sheet and rj is the magnetic diffusivity. If one consider 
the classical resistivity by Coulomb collision (Spitzer 1956), the Magnetic Reynolds number is 
enormously large (~ 10 14 ) in the solar corona. Therefore the reconnection rate is too small to 
explain the rapid energy release of solar flares, where the energy release occurs in 10-100 ta 
(= L/V A ; Alfven time). 

Petschek (1964) considered the effect of MHD slow shocks and proposed a model in which 
the reconnection rate is about 0.01-0.1 and almost independent of the Magnetic Reynolds number 
{Ma oc [lni? m ] -1 ). Several MHD simulations have shown that localized resistivity such as anoma- 
lous resistivity leads to Petschek type fast reconnection (e.g., Ugai 1986; Scholer 1989; Yokoyama 
et al. 1994). Erkaev, Semenov, & Jamitzky (2000) also showed that the Petschek regime was 
realized only when the resistivity is localized, by matching the Petschek-like solution in the inflow 
region and internal diffusion region solutions. Litvinenko, Forbes, & Priest (1996) examined the 
reconnection rate in flux-pile-up reconnection and suggested that the reconnection rate is slower 
than the value of Petschek reconnection by a factor of plasma beta. Nitta (2004) considered the 
self-similar solution of Petschek type reconnection (i.e., with slow shocks) in free space and found 
that the reconnection rate is ~ 0.05, nearly independent of the plasma beta and, possibly, of the 
magnetic Reynolds number, too. It is also worth noting that signature of slow mode MHD shocks 
associated with reconnection was found recently by X-ray observation of a giant arcade formation 
event (Shiota et al. 2003). 
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Observational measurement of reconnection rate is not straightforward because measurement 
of coronal magnetic field strength and inflow velocity is difficult. The only two direct observations 
of reconnection inflow reported in the literature are those of Yokoyama et al. (2001) and Lin et 
al. (2005). Yokoyama et al. (2001) found in EUV images two clear patterns approaching toward 
the reconnection X-point with a velocity about 5 km s _1 . The estimated reconnection rate was 
0.001-0.03. Lin et al. (2005) also found a signature of flows toward reconnecting current sheet in 
a CME-related event and obtained the reconnection rate of 0.01-0.23. 

Observations of the chromosphere provide us alternative ways to infer the reconnection rate in 
the corona. When magnetic reconnection occurs in the corona, the released energy is transported 
along the magnetic field to the chromosphere by nonthermal high energy particles and / or heat con- 
duction, brightening the footpoints of the reconnected field lines. Therefore successive brightenings 
in the chromosphere maps the successive magnetic reconnection in the corona. In the standard 
two-dimensional geometry, the successive reconnection of outer field lines results in the separation 
of two ribbons (Forbes & Lin 2000; Lin 2004), and the (dimensional) reconnection rate, that is 
also equal to the electric field in the reconnecting current sheet, is given by 

VinB cr = E cs = Vf oot Bf oot . (3) 

Here B cr and -B/ DO t are the magnetic field strength in the corona and at the footpoint, respectively, 
v foot is the separation velocity of the footpoints (two ribbons), E cs is the electric field in the current 
sheet, and Vi n B cr is the dimensional reconnection rate. Several authors derived E cs (or dimensional 
reconnection rate) using the observation of chromospheric two ribbons and photospheric magne- 
tograms (e.g., Qiu et al. 2002, 2004, Wang et al. 2003, Fletcher, Pollock & Potts 2003, Asai et al. 
2004b). 

Isobe et al. (2002) considered equation (3) and the energy release rate to derive the inflow 
velocity, magnetic field strength in the corona, and hence non-dimensional reconnection rate. They 
used soft X-ray images taken by the Soft X-ray Telescope (SXT; Tsuneta et al. 1991) aboard 
Yohkoh satellite to derive the energy release rate and Vf OQ t, and obtained Ma ~ 0.001 — 0.01 for 
the decay phase of a long duration flare. 

In this paper we apply the method presented by Isobe et al. (2002) to the impulsive phase of 
two ribbon flares and determine the reconnection rate Ma- We analyze three flares with different 
X-ray flux (one X-class, one M-class, and one C-class) to see whether there is a dependence of Ma 
on the flare class. We selected the flares that exhibit relatively symmetric and regular separation of 
the flare ribbons, because equation (3) is valid only in the approximately two-dimensional geometry 
such as shown in figure 1. We should note that it is possible that the successive brightening of the 
chromosphere is due to the successive magnetic reconnection along the magnetic neutral line, i.e., y 
direction in figure 1. In this case the velocity of the motion of chromospheric brightening cannot be 
used to calculate the electric field in the corona. The high electric field (90 V cm -1 ) in a complicated 
C9.0 class flare reported by Qiu et al. (2002) may be due to this effect. Indeed systematic successive 
reconnection along the neutral line has been observed in giant arcade formation events, where the 
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formation of X-ray loop progress along the neutral line with an apparent velocity of the order of 
10 km s" 1 (Isobe, Shibata, & Machida 2002). Interestingly, similar but faster (50 - 150 km s _1 ) 
evolution of reconnection along neutral line has been reported recently in a M-class flare observed 
by RHESSI (Grigis & Benz 2005). 

In order to calculate the energy release rate from SXT data as accurately as possible, we 
have carried out numerical simulations of a flare loop. We synthesized the count rates of different 
filters of SXT from the simulation result and calculate the thermal energy contents E t h that will 
be observed by the SXT, thus enabling us to convert the time derivative of E t h obtained from SXT 
data to actual energy release rate. 

In section 2, the method to determine the reconnection rate is briefly explained. In section 3, 
we describe the model and results of numerical simulations. In section 4 we describe the observation, 
data analyses, and results. Discussion and conclusions are given in section 5. 



2. HOW TO DETERMINE THE RECONNECTION RATE 

In this section we briefly summarize the procedure to calculate the reconnection rate, which is 
basically the same as that in Isobe et al. (2002) but uses different data set. 

We consider the energy release rate H by magnetic reconnection, that is given by the Poynting 
flux into the diffusion region; 

B 2 

H = 2— — Vi n Ly ;CS L Z:CS , (4) 

47T 

where L 2 CS is the vertical length of the reconnecting current sheet, and L ycs is the length of 
reconnecting current sheet along the magnetic neutral line. The geometry is illustrated in figure 1. 
The basic idea is that in equations (3) and (4) the unknown parameters are B cr and Vi n and the 
other parameters are obtained from observational data. 

As mentioned in the previous section, Vf OQ t can be measured from a time sequence of chro- 
mospheric images of two ribbons, and Bf oot can be measured from photospheric magnetograms. 
In this paper we use 1600 A images from the Transition Region and Coronal Explorer (TRACE; 
Handy et al. 1999) to measure vj 00 t, and magnetograms from Michelson Doppler Imager (MDI; 
Scherrer 1995) aboard the Solar and Heliospheric Observatory (SOHO; Domingo, Flech, & Poland 
1995) to measure Bf oot . 

Since the reconnecting current sheet is not usually visible, L y>cs and L z ^ cs are uncertain. We 
assume L ytCS is equal to the length of flare arcade L y . For L ZtCS , it seems reasonable to assume that 
in the impulsive phase L Z)CS is equals to the height of flare arcade L z . This assumption is consistent 
with MHD simulations (Chen & Shibata 2000; Shiota et al. 2003) and the inflow observation 
(Yokoyama et al. 2001). In the decay phase L z>cs may be much larger than L z . In this paper we 
analyze flares near the disk center, therefore we make further assumption that L z is equal to the 
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distance of the footpoints (flare ribbons) L x . It is not a bad assumption since many flare loop near 
the solar limb shows such geometry. Namely, we assume L 2jCS = L z = L x . 

To measure the energy release rate H, we use the SXT data. The temperature T and volume 
emission measure ey = n 2 V of the flare plasma can be derived from a pair of SXT images with 
different filters (Hara et al. 1992). Then density n and the thermal energy E^ can be calculated 
by assuming a suitable line-of-sight depth. The energy release rate H can be inferred from the 
time derivative of the thermal energy dE t h/dt, although the effect of radiative cooling, conductive 
cooling and chromospheric evaporation must be considered. In the next section we present the 
results of numerical simulations to convert the observed dE^jdt to the actual energy release rate 
H. 

Finally we calculate the unknown parameters B cr and Vi n from equations (3) and (4). The 
Alfven velocity va = B cr / \jAirm p nQ (m p : proton mass) is calculated using the calculated value 
of B cr and density no in the pre-flare phase derived from the SXT images. Thus we obtain the 
non-dimensional reconnection rate Ma = Vin/vA- To summarize, B cr ,Vi n , and Ma are given by 

Bcr = 2 lE , (5 ) 

v footBf 00 tL x Ly 

_ (vfootBf 00 t) 2 L x L y 

Vin- 2itH , W 

and 

Vin _ (VfootBf otf( L x L y ) 2 

v A ~ (2irH) 
where the parameters in the right hand sides are observable 



Ma = ^= v 1 V^h^o, (7) 



3. NUMERICAL SIMULATION OF A FLARE LOOP 

In order to derive the energy release rate H from SXT observations, we have carried out a set 
of numerical simulations of a flare loop with different loop length and energy input. The effects of 
Spitzer type thermal conduction, optically thin radiative cooling, and chromospheric evaporation 
are included. From the results of numerical simulations we calculated the expected count rates for 
specific filters of the SXT, and then re-calculated the average temperature, emission measure, and 
hence thermal energy by using the filter ratio method (Hara et al. 1992), as we usually do when we 
analyze the real observational data. By doing this we can calculate the ratio a of dE t h / dt obtained 
from the SXT data to the actual energy release rate H for various parameters; a = {dE t y t / dt) / H . 
Then a is used to calculate H from SXT data in section 4. In this paper we consider only the 
combination of Al 12 /im (A112) filter and Be 119 um (Be) filter. This pair of AU2 and Be filters 
has been frequently used in flare observations. 
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3.1. Basic Assumptions and Equations 

As a simple model of a flare loop, we consider a one-dimensional magnetic loop with semi- 
circular shape and constant cross section. The basic equations are: 

| + - 0, (S, 

dpv d 2 . dP 

-dt + dS {pv) = -^ii-sF' (9) 

de d d / dT\ 

- + -[(e + P)v] = -fa-)- mv - R + h , (10) 

where 

P = nk B T, e=-pv 2 + -. (11) 

z 7 — 1 

Here s is the distance along the loop from the photosphere. 7 = 5/3 is the ratio of the specific 
heats. g»(s) = gocos(irs/2L) is the acceleration by the gravity along the loop; L is the distance 
from photosphere to the top of the loop. The radiative loss function R = n 2 Q(T) is chosen to 
approximate the form for a optically thin plasma with normal solar abundances; Q(T) is shown in 
figure 2. The thermal conductivity is taken to be of the classical Spitzer form (Spitzer 1956), that 
is, K|| = k T 5 / 2 erg s 1 cm 1 K 1 and kq = 10 6 in cgs units. The resulting heat flow is assumed 
not to be flux-limited. The heating term h is divided into the static heating h s (s) and the flare 
heating hf(s,t). The static heating term h s (s) is assumed to balance the radiative losses in the 
initial condition at each point. The second term hf(s,t) is the flare heating, which is described 
later. The other parameters have their usual meanings. 

The basic equations are numerically solved as an initial-boundary value problem. The hydro- 
dynamic parts are solved by modified Lax-Wendroff method (Rubin and Burkstein 1967) and the 
heat conduction part is solved by the successive overrelaxation method (e.g., Hirsch 1988). The 
code has been extensively used for the modeling of solar flares (Hori et al. 1997), X-ray jets (Shi- 
mojo et al. 2001) and protostellar flares (Isobe et al. 2003). The detail of the numerical algorithms 
is described in the Appendix of Hori et al. (1997). 



3.2. Initial and Boundary Conditions 

The initial condition is hydrostatic atmosphere consisting of isothermal photosphere/chromosphere 
and hot corona. The temperature distribution in the initial condition is given by 

T = T pho + T -~ Tpfto {tanh (^») + 1}, (12) 

where T pho = 10 4 K, T cor = 2 x 10 6 K, s tr = 2.8 x 10 8 cm, and w tr = 10 7 cm. The initial density 
and pressure distributions are calculated from the temperature and gravity distributions by solving 
hydrostatic equation. 
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Since we assume that the static heating balances only with the radiative loss, the initial 
condition is not in strict energy equilibrium. Therefore, even without the flare heating slight 
decrease of temperature just above the transition region (s ~ st r ) and weak evaporation of the 
chromosphere spontaneously occur due to the conductive heat flux from the hot corona to the 
cold chromosphere. However, they have negligible effects on the result of simulation because the 
temperature increase and resultant evaporation by the flare are much stronger. The grid size As 
smoothly decreases from 3.2 x 10 6 cm at the photosphere to 2.0 x 10 7 cm around the transition 
region. In the corona, As increases from 2.0 x 10 6 cm to 1.0 x 10 7 cm. 

We assume symmetry about the loop top and solve only the half of the semi-circular loop, 
i.e., < s < L. Mirror boundary condition of the following form is imposed both at s = and 
s = L. Since our code requires that values of p, pv, and P are specified on one grid point at each 
boundary, we set pi = p2, P\V\ = P2^2, and Pi = P2 for the bottom boundary and pjy = Pn-i, 
Pnvn = Pn-iVn-i, and Pn = Pn-i for the top boundary. Here the subscripts denote the position 
of the grid (i.e., pi is the density on the ith grid point), and N is the total grid number. The 
bottom boundary condition does not affect the result of simulation because the gas and energy 
density near the bottom boundary is much larger than the upper atmosphere and therefore the 
bottom boundary is not perturbed significantly by the flare heating. We chose the mirror boundary 
because it is mathematically compatible and the total mass and energy in the simulation domain 
are conserved. 



3.3. Flare Heating 

The flare heating function hf(s,t) is of a form of spatially Gaussian, 

hf(s, t) = j^L exp[- (a " 2 y )2 ] (erg cm" 3 s" 1 ). (13) 

We adopt a = 2.0 x 10 s cm and Sfi r = L, i.e., the flare energy is deposited at the loop top. The 
total heat flux for the loop q(t) is set to be 

q (Q = qf{l-tanh{(t-tf)/w t f}) ^ ^ 

We adopt tf = 360 s and wtf = 1.7s, because the duration of the impulsive phase is about 6 
minutes for all the flares analyzed in this paper. We also performed simulations with uniform 
heating and footpoint heating to see how the result varies if different form of flare heating is used. 
The comparison of the different heating models is presented in the Appendix. 



3.4. Results of Numerical Simulations 



We performed the simulations with different loop length L and heat flux qf, L = 2.0, 3.0, 3.3, 
and 4.0 x 10 9 cm and qf = 0.4, 0.8, 4, 8, 40, 80, and 243 x 10 8 ergs" 1 cm -2 . The other parameters 
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were kept constant. In the following we show the result of the run with L = 3.3 x 10 9 cm and 
qj = 8.0 x 10 8 ergs -1 cm~ 2 . 

Figure 3 shows the temperature and density distributions at t = 0, 3.4, 34, and 340 s. Since 
the flare energy is injected at the loop top, the temperature at the loop top rises rapidly, and the 
energy is transported by heat conduction. The conduction front is seen near s = 1.2 x 10 9 cm 
at t = 3.4. When the conduction front reaches to the chromosphere, the cold and dense plasma 
is heated to the flare temperature, and the resultant high pressure drives the upward flows, i.e., 
chromospheric evaporation (t = 34 s). The evaporation plasma fills the flare loop and increases the 
soft X-ray emission (t = 340 s). 

From the temperature and density distributions in the result of simulation, the counts of signal 
that will be detected by SXT can be calculated by 

N 

^■ = E^( t *k 2As - ( 15 ) 
i=i 

where fj is the response function of the filter j, Ij is the count rate of SXT with filter j, i denotes the 
grid number, and Asj is the spacing of the ith grid. The left panel of figure 4 shows the synthesized 
SXT count rates for A112 filter {IaV- solid line) and for Be filter {I Be'- dashed line). Note that the 
unit is DN (data number) per second per unit area, where the area is that on the solar surface. 
Since one pixel of the full resolution image of SXT corresponds to 2.45 arcsec k, 1.8 x 10 8 cm on the 
solar surface, to obtain the count rate in DN s _1 pixel -1 , the values in figure 4 must be multiplied 
by 3.2 x 10 16 . 

From Iai and Ise, the average temperature T and emission measure €l = n 2 L in the flare loop 
can be calculated by the filter ratio method (Hara et al. 1992). Here is the emission measure of 
the loop with a cross section of 1 cm 2 . The thermal energy of the loop is given by 

e th = — — — r-£ = 3^1^ (erg cm" 2 ). (16) 
7-1 

The total thermal energy E t h is then given by E t h = e t hS, where S is the cross section of the flare 
loop(s) which is arbitrary in the simulation. Finally we obtain the ratio a; 

dEth/dt deth/dt 

« = 7} — = . ( 17 ) 

H q f 

where qf is the input heat flux of the flare. Note that the thermal energy calculated in this way is 
that of the high temperature plasma, because the SXT sees only the hot coronal plasma (T > 2 x 10 6 
K). The temporal variation of eth is shown in the right panel of figure 4. We defined the impulsive 
phase as 0.2tj/ r < 0.8t// r and calculated the average gradient of the ethif) in the impulsive phase by 
the least square fitting. The solid line in the right panel of figure 4 shows the result of least square 
fitting that yields deth/dt = 5.6 x 10 8 erg s _1 cm~ 2 . Since the input heat flux in this simulation is 
8.0 x 10 8 erg s _1 cm~ 2 , a = 0.7 in this case. 
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We calculated a for all the other parameters by the same procedure. The result is shown in 
figure 5. The left and right panels plot a against input heat flux and deth/dt, respectively. The 
numbers in the figure indicate the loop length. Although there are some exceptions, basically larger 
heat flux or shorter loop length result in smaller a. We interpret this tendency as follows. When 
the heat flux is larger or the loop length is shorter, the conductive flux into the chromosphere is 
larger, and hence the conductive cooling is more effective. Moreover, larger conductive flux leads to 
stronger chromospheric evaporation and larger density in the corona. Therefore radiative cooling 
is also effective in such cases. In the following analyses we use a to calculate the real energy release 
rate H from dE^/dt obtained from SXT images. 

4. RECONNECTION RATE 

Using the method described in section 2 and the result of the simulation presented in section 
3, we have analyzed three flares: an X2.3 class flare on 2000 November 24, an M3.7 class flare on 
2000 July 14, and a C8.9 class flare on 2000 November 16. These flares were selected because (1) 
they show relatively symmetric expansion of two ribbons in TRACE 1600 A image, (2) the SXT 
data taken with A112 and Be filters are available for the impulsive phases, (3) magnetograms from 
SOHO/MDI are available, and (4) the size (~ 4 x 10 9 cm) and the duration of the impulsive phase 
(~ 6 min) are similar. 

4.1. X2.3 Flare on 2000 November 24 

This flare is one of the three X-class homologous flares occurred on 2000 November 24 in NOAA 
9236. This active region produced many homologous flares and coronal mass ejections and has been 
extensively studied (e.g., Nitta & Hudson 2001; Zhang & Wang 2002; Zhang, Solanki, & Wang 
2003; Moon et al. 2003; Takasaki et al. 2004) A quantitative analysis of the ribbon separation 
and hard X-ray emission of the homologous flares was presented by Takasaki et al. (2004). 

The SXT image near the peak is shown in figure 6 (left panel). The solid lines indicate 
L x = 2.2 x 10 9 cm and L y = 3.5 x 10 9 cm measured by eyes. The effect of projection on the plain 
of sky is corrected. Since the X-ray arcade grows with time and the arcade and flare ribbons are 
not exactly symmetric, the measured L x and L y are rough estimates. In this paper we consider the 
spatial and temporal average of the reconnection rate, so we use above values of L x and L y as the 
representative values to calculate the reconnection rate. In order to study the spatial and temporal 
variation of the reconnection rate, the size and other parameters must be measured more carefully 
considering the complicated and asymmetric structure. 

The middle panel of figure 6 shows the temporal variation of the intensity integrated over the 
flare arcade. The solid and dotted lines are for A112 and Be filter images, respectively. Since the 
SXT instrument cannot obtain two images with different filters simultaneously, we calculate the 
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intensity of A112 filter images at the same time of Be filter images by linear interpolation, and 
then calculate the temperature T and the volume emission measure ey = n 2 V by the filter ratio 
method. We assume that the line of light length of the arcade (~ height of the arcade) is equals to 
L x , the distance between the footpoints. This assumption is justified because we empirically know 
from the observations of limb flares that the height of a flare arcade is approximately equals to the 
distance between its footpoints. We also assume the volume filling factor / = 0.1; V = L x L y f. 
This is also from our empirical knowledge that inner part of a flare arcade usually looks like a void. 
With these assumptions the thermal energy of the flare arcade is given by 

Eth = 3nk B TV = 3k B T^V. (18) 

The right panel of figure 6 shows the temporal variation of E t h- The solid line is the least square 
fitting in the impulsive phase. From the gradient of the solid line we obtain dEth/dt = 1.6 x 10 28 
erg s -1 . 

We use the result of numerical simulation to calculate the real energy release rate H from 
dE t h/dt. In order to compare with the one-dimensional numerical simulation, we need to know the 
energy release rate per unit area de t h/dt. To obtain this, we simply divide total energy release rate 
dEth/dt by the apparent area of the flare arcade L x L y ; deth/dt = (dE t h/dt)/(L x L y ) = 2.0 x 10 9 erg 
s _1 cm~ 2 . From the assumed geometry we estimate that the loop length (distance from footpoint 
to looptop) is about 1.3L X ~ 3 x 10 9 cm. Then from figure 5 we see that a = 0.6 for this flare. 
Thus we obtain H = (dEth/dt) /a = 2.7 x 10 28 erg s _1 . These parameters are summarized in table 
1. 

Next step is to measure the separation velocity Vf oot and field strength Bf oot of the flare ribbons. 
Although the three flares analyzed in this paper have relatively symmetric and simple ribbons, Vf oot 
and Bf 00 t are different at the different position on the ribbons. Since we are interested in the average 
values, Vfoot and Bf oot are measured by following procedure. (1) Two TRACE 1600 A images, one 
near the beginning of the impulsive phase and one just before the end of the impulsive phase, are 
selected. (2) The MDI magnetogram at the nearest time is selected and coaligned with the TRACE 
1600 A images. The coalignment is done by taking cross-correlation of the TRACE white-light 
image and the MDI continuum image. (3) The outer edges of the ribbons in the TRACE 1600 
A images are defined by visual inspection of the images (see figure 7). (4) The magnetic flux 
4> = J BdA, where A is the area swept by the outer edge of the flare ribbons, is calculated from 
the coaligned magnetogram. We assume that magnetic field is vertical at the photosphere, and 
correct the effect of projection. The average field strength Bf QO t is given by Bf QO t = <P/A. (5) The 
(average) separation velocity Vf QO t is calculated by 

A 

V foot = J 77 TTi (19) 

L •ribbon\ t 2 ~ t l) 

where L r ib& on is the length of the ribbon, and ti and t\ are the times of the TRACE images {ti > t\). 
This procedure is applied to each ribbon separately. We use the average Bf OQ t and v f QO t of the two 
ribbons to calculate the reconnection rate. 
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Figure 7 shows the TRACE 1600 A images near the beginning (t\ =15:07:34 UT) and end 
{p2 =15:14:21 UT) of the impulsive phase, as well as the coaligned MDI magnetogram. The dashed 
lines indicate the outer edges of the ribbons. The magnetic flux <p are calculated by summing up 
the flux of all the pixels between the ribbon edges at t = t\ and t = t^. The solid lines beside the 
locations of the ribbons (dashed lines) indicate L ri }^ on for the east (left on the image) ribbon and 
the west (right on the image) ribbon. 

Table 2 shows 4>, L^^on, Bf oot and Vf oot for the east and west ribbons and their average values. 
Effect of the projection has been corrected. It should be noted that the values of the magnetic 
flux 4> of the two ribbons are not equal. Such flux unbalance has also been reported by Fletcher & 
Hudson (2001). Finally the B 

en Vim Ma, and E cs are calculated by equations (3), (5), (6), and (7) 
using the average values of Bf ao t and Vi n as well as H, L x ,L y , and no derived from the SXT data. 
The results are shown in table 1. We obtained B c — 41 G, Vi n — 1.3 x 10 cm s , Ma — 0.047, 
and E cs = 539 V m _1 . 



4.2. M3.7 Flare on 2000 July 14 

This flare occurred in NOAA 9077 at 13:45 UT on 2000 July 14, the same day as the famous 
2000 Bastille Day flare but in the different active region. The SXT image, light curves, and the 
temporal variation of E t h are shown in figure 8 in the same manner as figure 6. Since the solid line 
showing the least square fitting of E t h(t) is almost invisible, a line with the same slope (= dE^/dt) 
is also shown in the figure. 

Following the same procedure as the 2000 Nov 24 flare, we obtained L x = 1.5 x 10 9 cm, 
L y = 4.2 x 10 9 cm, and dE^/dt = 3.5 x 10 27 erg s _1 . The ratio a is taken from figure 5. The loop 
length is 1.3L X « 2 x 10 9 cm and the energy release rate per unit area is (deth/dt) ~ 5.6 x 10 8 erg 
s _1 cm~ 2 , hence a = 0.55 and H = (dE th /dt)/a = 6.4 x 10 28 erg s _1 . 

The TRACE 1600 A images and magnetogram of the flare are shown in figure 9, and the values 
of 4>, Bf oot , and Vf oot are shown in table 3. When we define the outer edge of the flare ribbons to 
calculate the magnetic flux, the upper part of the east (left) ribbon and the lower part of the west 
(right) ribbon that exhibit irregular shape are neglected because they seem to deviate from the 
standard two-dimensional geometry. Nevertheless the unbalance of <j) is small. Finally we obtain 
B cr = 44 G, v in = 3.2 x 10 6 cm s -1 , M A = 0.015, and E cs = 143 V m" 1 . Other parameters arc 
shown in table 1. 



4.3. C8.9 Flare on 2000 November 16 

This flare occurred at 00:20 UT on 2000 November 16 in NOAA 9231. The SXT image, light 
curves, and the temporal variation of E t h are shown in figure 10, and the TRACE 1600 A images 
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and the MDI magnetogram are shown in figure 11, in the same manner as previous two flares. From 
the SXT data we obtain L x = 3.1 x 10 9 cm, L y = 5.2 x 10 9 cm, and dE th /dt = 1.2 x 10 27 erg s^ 1 . 
The loop length 1.3L X ~ 4.0 x 10 9 cm and the energy release rate per unit area deth/dt = 7.4 x 10 7 
erg cm -2 , hence from figure 5 a = 0.65 and H = 1.9 x 10 27 erg s _1 . Measurement of <p,Bf oa t, 
and Vf oot are also done in the same manner. The results are shown in table 4. Finally we obtain 
B cr = 11 (G), v in = 6.7 x 10 6 cm s _1 , M A = 0.071, and E cs = 70 V m -1 . Other parameters are 
shown in table 1. 



4.4. Estimate of Uncertainty 

It is difficult to estimate the uncertainties of the obtained values of MA,V{ n , and B cr , because 
there are many factors that cause the uncertainty, though our method includes fewer assumptions 
than previous studies. Here we try to make a rough estimation of the uncertainty. First of all we 
should note that the obtained reconnection rate is the spatial and temporal average in the impulsive 
phase. Actually the reconnection process is often noisy and intermittent, which is clear from the 
highly variable, bursty light curves of nonthermal emissions such as hard X-ray and microwave. 
Careful examinations of the flare ribbons also suggest that the reconnection rate is noisy both in 
space and in time (Saba, Gaeng, h, Tarbell 2002; Fletcher, Pollock, & Potts 2003). Therefore 
each elementary process of reconnection in the flares can be faster than the average value obtained 
here. 

The thermal energy contents of the flare calculated in this study is consistent with previous 
works (e.g., Ohyama & Shibata 1997, 1998; Emslie et al. 2004). However, we wish to measure the 
energy release rate H as accurately as possible, because uncertainty in H has significant influence 
on the calculated value of the reconnection rate. The possible source of the uncertainty in the 
energy release rate H is: (1) temperature and emission measure diagnostics by SXT, (2) error in 
the conversion from dEth/dt (measured from SXT data) to H using the results of the numerical 
simulations, and (3) measurement of the volume V. As for (1), the statistical error due to the 
photon noise is negligible. As for (2), a may depends on the detail of the simulation model such 
as the initial condition, location of the flare heating, etc. We performed simulations with different 
parameters and found that the detail of initial condition does not have significant consequence in 
the result when the total energy of the flare is much larger than the initial thermal energy contents 
of the loop. On the other hand, simulations with different heating locations indicate that if most of 
the flare energy is deposited deep in the chromosphere, a parameter may be significantly smaller. 
However, within the range of reasonable parameters the uncertainty in the ratio a is about ±20%. 
The result of simulations with different heating locations is presented in the Appendix. 

Probably the largest uncertainty in H comes from (3), the uncertainty of the volume. We 
made assumptions that the line of sight length of an arcade is equal to its footpoint distance L x 
and a volume filling factor of 0.1. We cannot tell precisely how good these assumptions are, but 
we empirically know that the assumption of the geometry (i.e., L z = L x ) is not a bad assumption, 
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and probably the uncertainty is a factor of 2 at most. An alternative way to estimate the line of 
sight length may be to use the scaling law that relates the loop length (from apex to footpoint) 
to the rise time, decay time, and temperature of the flare (Metcalf & Fisher 1996). We applied 
this scaling law by Metcalf & Fisher (1996) to the flares analyzed in this paper and obtained the 
loop length of 1.3 x 10 9 cm, 1.3 x 10 9 cm, and 2.6 x 10 9 cm for the X2.3, M3.7, and C8.9 flares, 
respectively. These values are 1.5-2.2 times smaller that our assumption that the loop length is 
equal to 1.3L X . This also supports that the uncertainty of the line of sight length is about a factor 
of 2. 

The uncertainty of the volume filling factor is more difficult to address, but it seems that 
an uncertainty of a factor of 5 is reasonable. The upper limit comes again from our empirical 
knowledge of the geometry of flare loops/arcades, and the lower limit comes from that if the filling 
factor is so small as 0.01, the number density of the plasma in the loops must be unreasonably 
large (~ 10 13 cm -3 ). 

Combining the uncertainties of the line of sight length and the filling factor, we conclude that 
the uncertainty of the volume is a factor of 10. From equation (18) we see that the thermal energy 
is proportional to VV, hence the uncertainty in E t h is about a factor of \/l0. Combining (l)-(3), 
we estimate that uncertainty in H is about a factor of 4. 

About the measurement of Bf oot and Vf oot , There must be errors that come from the alignment 
of the TRACE 1600 A ribbons and MDI magnetogram. The alignment of the TRACE images and 
MDI magnetograms has been done by taking the cross correlation of the white light images of 
TRACE and the continuum images of MDI. We found that the correlation between the two images 
was good in all the flares and estimated the error is about 1 MDI pixel. This is consistent with the 
work by Fletcher & Hudson (2001) who used the same method for the coalignment of TRACE 1600 
A images and MDI magnetogram. The magnetic flux <fi measured by shifting the magnetogram by 
1 pixel yield the uncertainty of a few %. 

However, the unbalance of <p between the two ribbons indicates that there must be other source 
of uncertainty. As discussed by Fletcher & Hudson (2001) in detail, the uncertainty may come 
from: (1) the inclination of the photospheric field (2) not all the flux below the illuminated ribbons 
are involved in the reconnection event (3) uncertainty of the sensitivity of MDI magnetogram, 
especially to the small scale, and so on. Although we do not know which factor is dominant, for the 
present purpose, it is enough to roughly estimate the uncertainty of the value of product Bf oot Vf 00 t. 
In the three flares analyzed in this paper, the largest difference of E cs = Bf oot Vf 00 t between two 
ribbons is about a factor of 2.4 (2000 Nov 24 flare). If we adopt the values for the east and west 
ribbons as the upper and lower limits, the uncertainty of Bf oot Vf 00 t may be considered to be a 
factor of y/2A rj 1.5; 



togw( B footVfoot) = log 10 (B f oot V f 00 t) average ± log 10 1.5 (« 0.2). 



(20) 
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From equations (5) - (7), we see 

B cr OC ^(BfootVfoot) 1 
V in OC H^iBfootVfootf 

M A oc H- 2 (B foot v foot )- 3 . 

Therefore, the estimated uncertainties are ± log 10 (4 1 x 1.5 1 ) ~ 0.8 for log 10 B cr , ± log 10 (4 1 x 1.5 2 ) ~ 
1.0 for log 10 v in , and ±log 10 (4 2 x 1.5 3 ) « 1.7 for log 10 MA- 

For the uncertainty of the reconnection rate Ma, the uncertainty of the prefiare density no 
should be also taken into account. The largest uncertainty of n$ comes from the uncertainty of the 
line of sight length I. However, since va oc n Q oc Z 1 / 4 , the uncertainty of the prefiare density can 
be neglected. 

5. DISCUSSION 

By considering the separation velocity of the flare ribbons and energy release rate, we have 
obtained dimensionless reconnection rate Ma as well as the inflow speed Vi n and magnetic field 
strength B cr in the corona for the impulsive phases of three flares. The obtained values of recon- 
nection rate (0.015-0.071) are consistent with the fast reconnection models (e.g., Petschek 1964). 
In other words we have confirmed that the motion of flare ribbons and the energetics of the flares 
are consistent with the standard scenario that the magnetic energy is released via fast magnetic 
reconnection. Furthermore, it seems that the large energy release rate H and electric field E cs in 
the X-class flare are simply due to the large magnetic field strength and Alfven velocity in the 
corona. 

Previously several authors tried to estimate the non-dimensional reconnection rate from ob- 
servational data (e.g., Dere 1996; Tsuneta 1996; Ohyama & Shibata 1997, 1998; Forbes & Lin 
2000) and obtained similar values. Our method has advantages to previous studies that it requires 
fewer assumptions and that the inflow velocity and magnetic field strength in the corona are also 
obtained (see discussion in Isobe et al. 2002). Therefore this method is also useful to estimate the 
coronal magnetic field strength. 

As mentioned in section 4.4, the reconnection rate obtained above is the spatial and temporal 
average in the impulsive phase. Observations have shown that actually reconnection in solar flares 
is quite intermittent, both in time and in space (e.g., Karlicky et al. 2005; Fletcher, Pollock, & 
Potts 2003). Numerical simulations have also demonstrated that reconnection is intrinsically time 
dependent (e.g., Tanuma et al. 2001). Furthermore, based on the three-dimensional simulation 
Isobe et al. (2005) suggested that reconnection may be (intrinsically) intermittent not only in time 
but also in space. Therefore, it is likely that in each elemental reconnection the reconnection rate 
is larger than the spatial and temporal average. In principle, the spatial and temporal variation of 
the reconnection rate can be determined by the same method used in this paper. In order to discuss 
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the spatial and temporal variation of the non-dimensional reconnection rate, accurate alignment 
of the soft X-ray images (to derive the energy release rate), chromospheric images (such as Ha or 
TRACE 1600 A), and the magnetogram are necessary. Although such analysis is possible using 
the present data, we hope that the data from the X-Ray Telescope (XRT) aboard Solar-B will be 
very suitable because of the high resolution and better temperature diagnostics. 

We obtained the inflow velocity Vi n = 30 — 130 km s _1 . Yokoyama et al. (2001) measured the 
velocity of the apparent motion of the bright structures in the EUV image of the inflow region and 
obtained Vi n ~ 5 km s _1 . This value is smaller than the those obtained in this paper, probably 
because Yokoyama et al. (2001) measured Vi n in the decay phase. Recently several studies related 
to Yokoyama et al. (2001) have been done. Narukage & Shibata (2005) found more examples of 
similar inflow events in SOHO/EIT data and obtained similar values of the reconnection rate. Chen 
et al. (2004) calculated the EUV image (FeXII 195 A emission line) from the result of numerical 
simulation and suggested the actual inflow velocity is larger than the velocity of the motion of 
bright pattern in the EUV image. On the other hand, Noglik, Walsh, Ireland (2005) analyzed the 
same event as Yokoyama et al. (2001) using basically the same method as this paper and obtained 
the reconnection rate that is consistent with Yokoyama et al. (2001). More accurate measurement 
of Vi n should be done by detecting the Doppler shift by spectroscopy. It will be an important target 
of the EUV Imaging Spectrometer (EIS) of the Solar-B, the solar observing satellite to be launched 
in 2006. 

In order to make a constraint on the theoretical models of reconnection, it is valuable to 
examine the relation between the reconnection rate and the magnetic Reynolds number. In the 
laboratory experiment, Ji et al. (1998) found that the measured reconnection rate in the 2D 
reconnection experiment can be explained by a generalized Sweet-Parker model that incorporates 
the compressibility and the effective resistivity. However, the magnetic Reynolds number of the 
reconnection experiments is less than 10 3 , hence it is not clear if the same model can be applied to 
the space and astrophysical reconnection where the magnetic Reynolds number is many orders of 
magnitude larger. 

The magnetic Reynolds numbers for the analyzed three flares are calculated using the vertical 
length of the current sheet L 2jCS = L x and the classical magnetic diffusivity given by (Priest 1982) 

7] = 5.2x 10 11 In AT -3/2 cm 2 s~\ (21) 

where for the Coulomb logarithm In A ~ 20 for the parameters considered here (T ~ 10 7 K, n ~ 10 9 
cm -3 ). To calculate r\ we use the temperature derived from SXT data outside the flaring arcade 
(To), as was done to derive no- 

The values of To and R m are shown in table 1. For the three flares analyzed here, there seems 
to be no dependence of the reconnection rate on the Magnetic Reynolds number. Of course the 
accuracy and the number of data is not enough to derive a conclusion on this issue. A statistical 
study of the reconnection rate of a large number of events will give us further insight into the 
physics to determine the reconnection rate. 
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A. Effect of Flare Heating Model 

The form of the flare heating used in this paper assumes that all the released energy is injected 
at the loop top as the thermal energy and transported to the chromosphere by thermal conduction. 
However, if the significant fraction of the magnetic energy is converted to the non-thermal high 
energy particles during reconnection, the released energy is transported by the high energy particles 
and therefore the heating of the flare loop occurs in the transition region or chromosphere where 
the kinetic energy of the particles is thermalized through collision. Possibly, the factor a derived 
from the simulations may differ if the location of the flare heating is different. In order to test 
how the result depends on the location of the heating, we have performed simulations with uniform 
heating and footpoint heating. 

The uniform heating function is given by 

hf{S}t) = *M ^H(s-Ztr)/w tr } + ^ (A1) 

? 2 

where £ = J L tanh[(s — zt r )/wtr]ds is the normalization factor. The footpoint heating function is 
given by 

hf(s, t) = eM JS ~l tr 2 +d)2 } (erg cm" 3 s" 1 ), (A2) 

where d is the depth of the energy deposition from the transition region, and <tj = 6x 10 7 cm. We 
examined the cases with d = 0, 2, 4, and 6 x 10 7 cm. Compared with equation (13), the factor of 
1/2 in equation (A2) is introduced so that the total heat flux in the loop is the same for the same 
value of q(t). 

Figure 12 shows the temporal variation of e^, the total thermal energy in the loop derived 
from SXT filter ratio method. The heat flux qf is fixed to 8 x 10 s erg s" 1 cm~ 2 in all the cases. The 
evolution of the uniform heating case is almost identical to the loop top case which is described in 
section 3.4. This is because the temperature is almost kept constant in the corona due to the very 
effective thermal conduction. Therefore spatial distribution of the flare heating has little effect as 
long as it is in the corona. 

On the other hand, footpoint heating models result in smaller observed thermal energy. As 
shown in figure 12, e t h (and hence de t h/dt) is smaller when d is larger. By examining the density 
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and temperature profiles we found that this tendency is because the larger density in the location 
of flare heating results in larger radiative cooling. In the case of d = and d = 2 x 10 7 cm the 
evolution of the flare loop is similar to that of loop top heating model, whereas in the case of 
d = 6 x 10 7 cm almost all the flare energy is immediately radiated in the chromosphere. The latter 
is probably an extreme case and not realistic, because in this case we expect strong hard X-ray 
emission by the precipitation of high energy particles and very weak soft X-ray emission, that is not 
consistent with the observed correlation with hard X-ray and soft X-ray intensities. Furthermore, 
by carefully examining figure 12 we find that difference in deth/dt (gradient of the curves in the 
rising phase) is smaller than that of the peak values of e^, except for the d = 6 x 10 7 cm model. 
Therefore we conclude that the factor a derived in this paper assuming the loop top heating is 
quantitatively reasonable, with uncertainty of about ±20%. 

We should note that the large radiative cooling in the chromosphere may be due to the opti- 
cally thin radiative cooling model, which is probaly not a good approximation in the actual solar 
chromosphere. More realistic modeling with proper treatment of non-LTE radiative transfer (e.g., 
Abbett & Hawley 1999) is necessary for more acurate evaluation of the energy release rate. Fur- 
thermore, careful modeling of the chromospheric and transion region emissions is desirable to make 
spectroscopic diagnostics to determine what fraction of the released energy goes to the non-thermal 
particles and what goes to the thermal energy. 
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Fig. 2. — Radiative loss function Q(T). 
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Fig. 3. — Temperature and density distribution of the simulation result with typical parameters 
(see the text). The labels 1, 2, 3, and 4 in the figure indicate the time; t = 0, 3.4, 34, and 340 s, 
respectively. 




Fig. 4. — Left: Synthesized light curves. The solid and dashed lines are the light curves of Thick 
Al. and Be filters, respectively. They must be multiplied by 3.2 x 10 16 to obtain the value in DN 
s _1 pixel -1 for full resolution image. Right: Temporal variation of e^, the thermal energy per 
unit area derived by the SXT filter ratio method. The solid line is the least square fitting of the 
impulsive phase. 
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Fig. 5. — The ratio a = (de t h/dt)/qf plotted against qf (left) and de t h/dt (right). The numbers in 
the figure indicate the loop length (fO 9 cm). 
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Fig. 6.— Left: SXT image of the X2.3 flare on 2000 Nov. 24. The image is negative. The solid 
lines indicate L x (shorter line) and L y (longer line) measured by eyes. Middle: Light curves of the 
Al 12 filter images (solid line) and Be filter images (dotted line). Right: Temporal variation of 
thermal energy E^. The solid line indicates the least square fitting of the impulsive phase, that 
yields dE th /dt = 1.6 x 10 28 erg s -1 . 
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Fig. 7.— Left and middle: TRACE 1600 A images of the X2.3 flare at h =15:07:34 and 
t2 =15:14:21. The dashed lines indicate the outer edge of the flare ribbons at each time. Right: 
MDI magnetogram of the same field of view. The dashed lines indicate the locations of the outer 
edges of the ribbons both at t = t\ and t<i- The solid lines in the left panel indicate the length of 
the ribbon L r ibbon f° r the east (left) and west (right) ribbons. 
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Fig. 8.— The same figure as figure 6 for the M3.7 flare on 2000 Jul. 14. dE th /dt = 3.5 x 10 27 erg 
s- 1 . 




Fig. 9.— 



The same figure as figure 7 for the M3.7 flare on 2000 Jul. 14. 
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Fig. 11. — The same figure as figure 7 for the C8.9 flare on 2000 Nov. 16. 



Table 1: Parameters of the flares. 



Parameter 


2000 Nov 24 


2000 Jul 14 


2000 Nov 16 


GOES class 


X2.3 


M3.7 


C8.9 


H (erg s _1 ) 


2.7e28 


6.4e27 


1.9e27 


v foot (cm s" 1 ) 


1.2e6 


1.2e6 


6.7e5 


Bfoot{G) 


449 


117 


106 


L x (cm) 


2.2e9 


1.5e9 


3.1e9 


L y (cm) 


3.5e9 


4.2e9 


5.2e9 


no (cm -3 ) 


1.0e9 


2.0e9 


6.e8 


7b (K) 


1.0e7 


8.0e6 


8.e6 


B C r (G) 


41 


44 


11 


v in (cm s _1 ) 


1.3e7 


3.2e6 


6.7e6 


va (cm s _1 ) 


2.8e8 


2.1e8 


9.4e7 


M A 


0.047 


0.015 


0.071 


Rrn 


1.8el5 


7.0el4 


6.3el4 


E cs (V m- 1 ) 


539 


143 


70 


a 


0.6 


0.55 


0.65 
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Fig. 12. — Temporal variation of eth in the simulations with various flare heating fucntion. Plotted 
are the loop top heating, the uniform heating and the footpoint heating with d = 0, 2, 4 and 6 x 10 7 
cm, which are indicated in the figure. 

Table 2: Magnetic flux (ft, ribbon length L r ibbons, magnetic field strength Bf QO t and separation 
velocity Vf OQ t of the 2000 Nov. 24 flare. The effect of projection was corrected. 
Parameter east ribbon west ribbon average 
(ft (Mx) 1.0e21 7.6e20 8.8e20 

L ribbon (cm) 3.1e9 5.8e9 4.5e9 

B foot (G) 574 324 449 

Vf 00 t (cm s" 1 ) 1.4e6 1.0e6 1.2e6 



Table 3: Magnetic flux (ft, ribbon length L r ibbonsi magnetic field strength Bf oot and separation 

velocity Vf oot of the 2000 Jul. 14 flare. The effect of projection was corrected, 

parameter east ribbon west ribbon average 

(ft (Mx) 1.6e20 1.5e20 1.6e20 

L ribbon (cm) 3.4e9 6.6e9 5.0e9 

B foot (G) 144 91 117 

v foo t (cm s" 1 ) 1.4e6 1.0e6 1.2e6 
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Table 4: Magnetic flux (ft, ribbon length Ivi(,& ons , magnetic field strength Bf oot and separation 
velocity Vf OQ t of the 2000 Nov. 16 flare. The effect of projection was corrected. 



Parameter 


east ribbon 


west ribbon 


average 


(ft (Mx) 


1.4e20 


2.2e20 


1.8e20 


^ribbon ( cm ) 


5.3e9 


7.6e9 


6.4e9 


B foo t (G) 


143 


69 


106 


Vfoot (cm s" 1 ) 


4.0e5 


9.3e5 


6.7e5 



